arXiv:hep-ph/0211450vl 28 Nov 2002 


Pion Structure Function and Violation of 
the Momentum Sum Rule 

J.P. Lansberg* F. Bissey, J.R. Cudell, J. Cugnon, M. Jaminon and 

P. St assart 


Universite de Liege, Departement de Physique B5, Sart Tilman,B-4000 LIEGE 1, Belgium 


Abstract. 

We present a method to evaluate the pion structure functions from a box diagram 
calculation. Pion and constituent quark fields are coupled through the simplest pseudoscalar 
coupling. The 7 * 7 r —> qq cross-section is evaluated and related to the structure functions. 
We then show that the introduction of non-perturbative effects, related to the pion size 
and preserving gauge invariance, provides us with a straighforward relation with the quark 
distribution. It is predicted that higher-twist terms become negligible for Q 2 larger than 
about 2 GeV 2 and that quarks in the pion have a momentum fraction smaller than in the 
proton. We enlarge the discussion concerning this violation of the momentum sum rule, 
emphasizing that the sum rule is recovered in the chiral limit and also when the finite size 
condition is not imposed. 


1. INTRODUCTION 

The usual way to describe quarks inside hadrons relies on the structure functions. 
Quantum Chromodynamics (QCD), being the theory of the interactions between 
quarks and gluons, is supposed to predict the structure functions, but in fact, these 
quantities result from non-perturbative effects and hence are not easily modeled. 
QCD is only able to predict their evolution as the virtuality of the probing photon 
( Q 2 ) changes. It is our purpose here to build a simple model for these quantities. 

Phenomenological quark models, which possess some non-perturbative aspects 
and which are rather successful in reproducing low-energy properties of hadrons, are 
expected to help us understand the connection between deep-inelastic scattering 
(DIS) data and non-pertubative inputs. For very simple systems such as pions 
and other low-mass mesons, there exist some effective models ( e.g. the Narnbu- 
Jona-Lasinio (NJL) model) that incorporate, in some simplified way, special QCD 
features such as spontaneous chiral symmetry breaking and anomalies. 

However, investigations of the structure functions along these lines have given 
rather different results [Sj. ||, [|]. This situation originates from the fact that the 
NJL model needs to be regularized and that different regularizations yield different 
results. 

We present here an alternative way to evaluate of the pion structure function 
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FIGURE 1. Feynman diagrams for the process 7*7ro —> qq. 


based on Ref. [| 1 | , where the qqn vertex is represented by the simplest pseudoscalar 
coupling, i.e. ig^/ 5 and the pion size is mimicked by the introduction of a gauge 
preserving cut-off. 


2. CROSS-SECTION CALCULATION FOR 7 *7r 0 > QQ 
BY INTEGRATION OVER T 

In this section, we calculate the cross-section for 7 * 7 r 0 —> qq which can be related to 
the form factors W\ and W 2 as we shall see later. Another way to achieve this goal 
consists in calculating the imaginary part of the forward amplitude for 7*7r 0 —> 7*717. 

At the leading order in the loop expansion, we have the two diagrams shown 
in Fig. |l| which, when squared, give two types of contributions. The first one 
(corresponding to box-diagrams for the 7*717 —> 7*717 amplitude) enables us to 
calculate the structure functions straightforwardly. The second involves crossed 
quark lines and is related to processes where the two photons (in the 7*717 —> 7*717 
scattering) are connected to different quark lines. This is incompatible with a 
probabilistic interpretation in terms of structure functions. As we will conclude 
from the calculations, contributions of this kind are fortunately higher-twists, i.e. 
they scale as ^2, when we introduce the cut-off. 

Their introduction is nevertheless mandatory, as this ensures gauge invariance 
(see Ref. [IJ). This is required for a proper evaluation of the cross-section for the 
7*717 —*► qq process, which is a physical process, observable and thus gauge invariant. 


2.1. Kinematics of the Process 

For the sake of convenience, we work with cartesian coordinates in the CM 
(center of momentum) frame. The four-momenta of the four particles are given by: 


q= ( E(=q 0 ),q) J p= {E^,-q),k 3 = {k 3fi ,k 3 ),k 4 = (fc 4 , 0 ,~h). 


( 1 ) 
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FIGURE 2. Kinematics and definition of the angle 9 in the CM frame. 


Let us now introduce the Mandelstam variables: 

s = (p + q ) 2 = (k 3 + Aq ) 2 = Eq M — (E + E n ) 2 , 
t= (fc 3 -q ) 2 = ( k 4 -p) 2 ,u = ( k 3 -p ) 2 = (. k 4 -q ) 2 . 

We then define the square momenta: 

2 sV2 2 2 1 2 2 1 2 2 

q = —Q ,p = k 3 —m ,k 4 — m 


( 2 ) 
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The variables of Eq. ([[]) can be written in terms of these quantities and of 
Mandelstam variables. One has: 


k3,o = h,o = ^Y , N = M = 


( 4 ) 
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The angle 6 between the vectors k 3 and q (see Fig. |zj) is related to t by 

+ 2,2 , r s ~ m l-Q 2 , O ft ~y/(s-ml-Q 2 ) 2 + 4 sQ 2 

t = m +rn^ — s + Vs - ^ - \- 2 J- — m z — --——-cosd 
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Eventually we obtain: 


t - m 2 -?f + §+<f + f + f 

COS0 = 2 2 . _ 2 
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The physical values of t are lying between boundaries, corresponding to cos# = 
± 1 : 
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Finally, assuming the incident direction along the 0 axis and the scattering in 
the y — z plane, we get: 


2 Vi 


9 = ( s -’:t: g2 .o.o.d (s ~ m -~ Q2)2+4sq V (g.o,o,ysHQ?), 
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p = Uj m 2 | {s-m 2 n -Q 2 ) 2 + ^sQ 2 Us-ml-Qiy + AsQ^ 
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h = (^A-J^-rn 2 8m(6),J^-m 2 co8d), 


h = (^y, 0 ,J^-m 2 sm{ 0 ),-J^-m 2 cos 6 ). 


( 9 ) 


These quantities will be explicitly used in the determination of scalar products 
involving the polarisation vectors. 


2.2. Amplitude Calculation 


Applying Feynman rules for the processes shown in Fig. |l], we obtain the 
analytical expression of the square amplitude in a straightforward way. For a given 
polarisation of the virtual photon, i — L, T \, T 2 , and with well-known conventions^ 
we have for the squared of the first diagram shown in Fig. [|, 

\Mi \ 2 = 3 g 2 (e 2 u + e 2 d )Tr + m)^ 

^^ + m ^ \k 3 3 -q) 2+ -Zi 2 ^) £i ^" (10) 


The expression for the square of the second diagram as well as for the interference 
terms have similar expressions. 

We first proceed to the calculation of the expressions of the scalar products of the 
polarisation vectors with q,p,k$ and k 4 for the three possible polarisation states. 
We shall use the following definition of the longitudinal polarisation vector, 






{q z , 0,0 ,e). 


( 11 ) 
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The isospin/charge factor (e^ + e^) corresponds to the following choice of the isospin matrix: 
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One has 


£L-q = 0 , 


£l-P = 






£l ■ h = 
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■m 2 ~^ + i + £ 
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2 2 s/W+Q 2 

For the linear transverse polarisation vectors, we use the following definition: 

£ Ti ,/2 = ( 0 , 1 , 0 , 0 ), ^=( 0 , 0 , 1 , 0 ). ( 13 ) 

This and Eq. (|9|) lead to: 

£ T .-q = 0 = e Ti -V , * = 1,2, 

£t x ■ k$ = 0, 


£t 2 ' ^3 — 


-m * x— i -^— , 

4 \ V 2 v / F^v / ^+g^ 
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(14) 


2.3. Polarised Cross-Sections 


The total cross-section is given by the following relation ||, with i = L, T \, T 2 : 


do* 

d£ 


167tA 


|M:| 


(15) 


with[] A = A(s, —Q 2 ,m 2 ). The latter directly links the amplitude -which is obtained, 
in our case, with REDUCE- and the cross-section differential upon t. Integrating 
over t from f min to f max (see Eq. (^)) gives the total cross-sections for given 
polarisations: 


°X = 3g 2 {e 2 u + e d) J dt (16) 

with T^' = Tr ((-^ 4 + m) 7 5 7^ W 3 + Ht" (£"^ 2 7 5 ) • 

It is also of interest to consider 

Vcontr. = Y^x 3 9 2 (e 2 u + e 2 d ) J d t T IHl '(g IHl , + ^f) (17) 


since it does not involve any polarisation vectors and its calculation is somewhat 
easier. However, it should not confused with cr L -|-cr Tl +&t 2 ( see later), which could 
be measured in some limiting experimental conditions. 


3 The function A {x,y,z) is defined as A(x, y , z) = x 2 + y 2 + z 2 — 2 xy — 2 xz — 2 yz 






















2.4. The Cut-Off 


2.4-1- Cut-off and Higher-twists 

We introduce a cut-off in our calculations: instead of integrating over t from 
fmin to f max , we integrate f'rornQ Max{t m i n ,A t } to t max . As shown in Ref. [|I|], this 
is equivalent to setting a cut-off for the square of the relative momentum of the 
quarks inside the pion, which can be considered as resulting from the finite size of 
the pion. 

Irrespectively of the factors arising in the extraction of the structure functions 
(see section ([|)), we can simply convince ourselves that in the presence of a cut¬ 
off the contribution of the crossed diagrams (the first photon hits a quark and 
the second an antiquark) are higher-twists, by plotting the ratio of these two 
contributions as a function of Q 2 . The suppression is clear on Fig. ^ for whatever 
(finite) value of the cut-off; the curves are even indistinguishable in this case. 

As we already said, in order to get a physical interpretation of our results -in 
terms of structure functions-, we need the crossed diagram contributions to be 
suppressed. The importance of a cut-off, which is not a regulator in the studied 
process, is thus manifest. 



Q 2 (GeV 2 ) 

FIGURE 3. Ratio of the direct contributions to the crossed contributions to the cross-section 
for 7*7T° —> qq at x = 0.1 and for three values of the cut-off A. 
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Recall that t < 0. 






2.4-2. Cut-off and Gauge Invariance 


In most cases, the presence of a cut-off implies the breakdown of gauge invariance. 
This is the main reason why one commonly uses other regularization methods such 
as Pauli-Villars or dimensional regularization. In the context of an effective theory, 
a simple momentum cut-off has the critical advantage that the regularization 
parameter (which is kept to a finite value) can be interpreted as a physical one. In 
the case of dimensional regularization, the signification of fractional dimensional 
space is not at all clear. As for Pauli-Villars regularization, the introduction of 
wrong-sign fermionic contribution is not intuitive. 

Nevertheless, in our model, as the intermediate states in 7*7r° —> 7*7T° are on- 
shell, the cut-off on the internal loop momentum^] acts on a physical observable, 
that is t of the process 7*7r° —> qq. The conclusion is thus staightforward: this cut-off 
does not break gauge invariance. 

This procedure would have been fully justified if we had only one pion vertex, 
and thus one condition to implement. But we have two vertices with two different 
expressions for the relative momenta. To keep gauge invariance, the cut-off on t 
must be on the overall cross-section as the physical observable is the full cross- 
section (the square of the sum of the two diagrams) at a given t. 

To keep the constraint coming from the pion size, we should have imposed two 
conditions on the cross-section to be integrated over t. However, this would have 
given zero. As a makeshift we imposed one condition or the other. 

Even if the physical interpretation is not so obvious, we thus kept gauge invari¬ 
ance, information from the pion size, suppressed the crossed diagrams and hence 
obtained an interpretation in terms of structure functions. 


3. Wi, W 2 , F\ AND F 2 


At this stage we have calculated the polarised cross-sections and we should be able 
to link them to the structure functions for neutral pion. Before doing this, we shall 
Erst recall some mandatory formulae derived from the analysis of deep inelastic 
proton scattering. 

We introduce the usual form of the hadronic tensor 


W^(p,q) = {-g^ 


Q 2 


p-q. 


f p-q. 


)W 1 {Q , u) + \p„ + (—)qf \p v + (— )q v ] 


{ w 2 (Q 2 




Q 2 


Q 2 


M 2 


,( 18 ) 


with M, the hadron mass, and Wi(Q 2 ,u) and W 2 (Q 2 ,v) the structure functions. 

In the same spirit as for the definition of the cross-section of —> X 0, we 

define the polarised cross-sections for 7*7r 0 —> X by [j: 


5 which corresponds, up to some constants, to a cut-off on the relative mometum between the 
two quarks at the n°qq vertex. 

6 This definition is independent of the hadron spin because, on the one hand, we have to sum 
the contributions coming from the different polarisations and on the other hand, we have a spin 
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Defining the transverse cross-section through 


we can easily show that: 
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( 21 ) 


These last two formulae establish the link between the functions W\ and W 2 and 
previously calculated polarised cross-sections for 7*7r 0 —> qq (Eq. fM 


factor arising in the relation between \V llv and the hadronic current. 

' Besides we would like to stress that one has to be cautious if one wants to consider non¬ 
polarised processes. Usually, a simple contraction of the amplitude for a given helicity with the 
expression, 


~9 


[IIS 


q^q v 


( 22 ) 


seems to be the most convenient way to proceed, instead of considering separately the different 
polarisations. This is in most cases licit because 


i=L,T lt T 2 


eP S " = -g^- 


IsP 


(23) 


for massive particles, and thus the -strict- sum over polarisations gives the correct tensor. 
For spacelike virtual particles, nevertheless, we have the following normalisations 


£ t -£t = — 1 , 

£l'£l = +1. (24) 


Therefore, the following relation -which indeed reduces to Eq. (23) for massive particles for 
which Eq. (gj) does not hold-, 
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leads, in this case, to 
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but not to ctl + 2<jt ! 


















With all these tools, we can now extract the functions W\ and W 2 • Inverting for 
W\ and W 2 , we get 
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Vx 


87r 2 m 7r a 


<7 t 


and 


W 2 = 


Va 


2 xml 
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We then get F\ and F 2 through the usual scaling relations: 

F 1 (x) = MW 1 (Q 2 ,v) ,F 2 (x) = ^-W 2 (Q 2 , l y) ,x = -^ - , 

M 2p-q 2v 


(27) 


(28) 


= £ (29) 


still with M for the hadron mass. 


4. SUM RULES, PLOTS AND RESULTS 
4.1. Number of Particles Sum Rule 

The structure functions can now be related to the (valence) quark distributions: 

4 1 - 

Fi = —(u v (x)+u v (x)) + —(d, v (x) + d v (x)), (30) 

lo lo 

F 2 — 2xF\. (31) 

We stress that the last relation, known as the Callan-Gross relation, comes out 

of our calculation. This does indicate that our approximations have been done 
consistently^ 

So far, we have only described the model. In order to make predictions, we need 
to fix its parameters, namely A, m and g. The latter can be thought of as the 
normalisation of the quark wave function, and we determine it by imposing that 
there are only two constituent quarks in the pion. In our model, the valence quark 
distributions are equal 

u v (x) = u v (x) = d v (x) = d v (x) = v(x). (32) 

The condition v(x)dx = 1/2 then yields 

U (x)dx — (33) 


8 For charged pions, the additional diagrams implying a direct coupling of the virtual photon to 
the pion are suppressed by a factor 1/s and the leading-twist results are the same, except for the 
charge coefficients entering Eq. (^) . 







FIGURE 4. Values of the coupling constant g{Q 2 , A) which fulfill the sum rule (Eq. ©), for- 
the values of the parameters indicated at the top. 


As F\ is a function, not only of m and A, but also of Q 2 , this gives us a coupling 
constant that evolves with Q 2 . The resulting values of g are shown in Fig. [| (right): 
for a finite cut-off A, the cross-section at fixed g would grow with energy, until the 
pion reaches its maximum allowed size, in which case the cross section would remain 
constant. If we impose relation Eq. d33|) , this means that g(Q 2 ) will first decrease 
until the cut-off makes it reach a plateau for Q 2 !§> A 2 (in practice, the plateau is 
reached around Q 2 ps 2A 2 ). The plateau value depends on m/A (and m^/A). 

4.2. Momentum Sum Rule 

To constrain further our parameters, we can try to use the momentum sum rule 

/*! jo r i 

2 hr) =4 / xv(x)dx = — / F 2 (x)dx. (34) 

Jo 5 Jo 

In the parton model, this integral should be equal to 1 as Q 2 —> oo, as we do not 
have gluons in the model. 

Analyzing Fig. ||, it is clear that - at least for a finite value of A- the momentum 
sum rule is violated. For instance the momentum carried by the quarks at large Q 2 
is 0.6 for m = 300 MeV and A = 750 MeV. 

Nevertheless, there are two limits of our model that fulfill the condition 2 ( x) = 1 
automatically and for which the quarks carry the entire momentum of the pion. 
First of all, if we do not impose that the pion has a finite size, then asymptotically 
we get 2 (x) =1. In the limit = 0, we obtain the simple expression 

41n(2 v/m 2 ) — 3 
41n(2 u/jti 2 ) — 1 


2 (x) (A = oo, = 0) 


(35) 













FIGURE 5. Momentum fraction of the quarks inside the neutral pion (Eq. 34), as a function 
of Q 2 , for the values of the parameters indicated at the top of the figure. 


This means that in that regime, for sufficiently high Q 2 , the quarks behave as free 
particles, and the usual derivation based on the OPE holds ||. 

The second case where this holds confirms this interpretation: if we impose 
m < m, r /2, the expressions we have given develop an infrared divergence, which 
corresponds to the case where both quarks emerging from the pion are on-shell and 
free. This divergence can be re-absorbed into the normalisation Eq. ( fT3| ) of g, and 
the sum rule 2{x) = 1 is again automatically satisfied at large Q 2 . 

However, in the physical pion case, it makes more sense to consider that one 
of the quarks remains off-shell: the introduction of a cut-off changes the sum rule 
value, as the fields can never be considered as free. Hence, because of the Goldstone 
nature of the pion, one expects that the momentum sum rule will take a smaller 
value than in the case of other hadrons. This may explain why fits that assume the 
same momentum fraction for valence quarks in protons and pions || do not seem 
to leave any room for sea quarks [10 . 

We come back to the results for 2(x) shown in Fig. Again, the curves show 
a plateau at sufficiently large Q 2 , with a value depending on m/A and m n /A. It 
is easy to show that this value is always smaller than 1 as it can be guessed from 

Fig. |. 

Following the discussion in J[J, we choose the following conservative values for the 
remaining parameters of our model: m = 300 MeV or 360 MeV and A ~ 800 MeV. 
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FIGURE 6. Asymptotic value (for large Q 2 ) of the momentum fraction of the quarks inside 
the neutral pion (Eq. [34]), as a function of the cut-off A. 

4.3. Plots of F 2 




X X 

FIGURE 7. Structure function F 2 for the neutral pion. On the left, Q 2 =2 GeV 2 and the 
values of the parameters m and A are as indicated. On the right, m=0.3 GeV and A=0.75 GeV 
and the values of Q 2 are as indicated. 

Let us finally examine the properties of the distribution v(x), or equiva¬ 
lently, of the function F 2 . Some of our results are summarized in Fig. [7], for 
Q 2 = 2 GeV 2 . The most striking feature is the vanishing of this function for x 
larger than some value x max - This is again due to the fact that, because the quarks 
are not free in this model, the actual value of their mass does matter. It can clearly 
be seen that this effect originates from kinematical cuts (see Ref. 1]). 















As a result, putting the cut quarks on-shell will be easier for small than for large 
x. Therefore, the ^-distribution (iq) is expected to be enhanced on the low x side, 
leading to a momentum fraction smaller than unity. The vanishing at large x is 
not obtained in similar works, in particular in the one of Ref. j2j. In this reference, 
the Bjorken limit is taken first and the kinematical constraint is not applied. This 
procedure seems questionable for evaluating cross-sections at finite Q 2 . 


5. CONCLUSION 

We have discussed the simplest model allowing to relate 7 * 7 r —> qq to the pion 
quark distributions. The model consists in calculating the simplest diagrams for a 
7 5 coupling between pion and constituent quark fields. 

We show that the crossed diagrams, which are to be included in order to 
guarantee gauge invariance, preclude the existence of the relation between the 
the cross-section and the quark distribution. The introduction of a cut-off on t 
allows such a relation without breaking gauge invariance. The cut-off on t being 
equivalent to a cut-off on the square of the relative momentum of the quarks inside 
pion, this effect can be traced back to the finite size of the pion. For reasonable 
value of the cut-off, the crossed diagrams become negligible as soon as Q 2 is larger 
than 2 GeV 2 allowing quark distributions to be meaningful even for relatively small 

Q 2 - 

One of the most important outcomes of our calculation is the deviation from 
the momentum sum rule. However, we recover this sum rule ( 2(x ) = 1) in two 
limiting cases, which correspond to free quarks. The first one is the limit A —> 00 . 
The second one is for m < Then, an infrared divergence appears, which can 
be reabsorbed into the normalisation of g and the momentum sum rule is again 
recovered for large Q 2 . Hence, the violation of the sum rule can tentatively be 
attributed to the Goldstone boson nature of the pion. 
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